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Q ■ We investigate the entanglement dynamics and coherence of two two-level atoms 

. interacting with two coherent fields of two spatially separated and dissipative cavities. 

It is in particular shown that entanglement sudden death is obtained within a short 

Qh' interacting time. However, after a long interaction time a long-lived entanglement 
+i> ■ 

j5 ' is shown, that is, the initial entanglement of two atoms could be partially preserved. 

^ , , In addition, the coherence of the two atoms will not be lost during the evolution. 
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I. INTRODUCTION 



- ^ , Quantum entanglement is at the heart of quantum information processing and quantum 

' computation 

mm 

. In recent years, many efforts have been devoted to the study of 



the evolution of joint systems formed by two qubits 



In particular, Yu and Eberly 
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have found out that the single-qubit dynamics and the 
global dynamics of an entangled two-qubit system subjected to independent environments 
are rather different. For a single-qubit system subjected to an environment, the local coher- 
ence decays asymptotically while the entanglement of an entangled two-qubit system may 
disappear for a finite time during the dynamics evolution. The nonsmooth finite-time dis- 
appearance of entanglement is called "entanglement sudden death" (ESD). More recently, 
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the ESD phenomenon has been observed in the laboratory by two groups [ly, [iTj. 

One the other hand, the long-hved entanglement in cavity QED or solid state systems was 
investigated by several authors 
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221 ] ■ In Ref. [18|, the authors have discussed 



the effects of the classical noise on entanglement dynamics and pointed out that there is 
long-lived entanglement in the presence of global dephasing noise. The influence of white 
noise on the entanglement dynamics of two atoms within two cavities has been investigated 
by Xu and Li |l^. They found out that if only one atom is driven by the white-noise field, 
then there is long-lived entanglement. The entanglement dynamics in two effective two- level 
trapped ions interacting with a laser field has been studied in [2^. It has been shown by 
Dajka, Mierzejewski, and Luczka that there is long-lived entanglement of some composite 
systems even if coupled with a thermal bath |2l|. It has also been pointed out that phase- 
damped cavity can lead to long-lived entanglement in a quantum system consisting of a 
single-Cooper-pair box irradiated by a quantized field |22 |. 

In the present paper, we investigate the entanglement dynamics and coherence of a quan- 
tum system formed by two two-level atoms within two spatially separated and dissipative 
cavities in the dispersive limit by the employing concurrence and linear entropy, respec- 
tively. The two atoms are initially entangled and the cavities are initially prepared in 
coherent states. We show that the ESD phenomenon appears in the present system and 
the entanglement of two atoms decreases with the time in the short-term. However, the 
long-term behavior is very different since a long survival of the entanglement is shown in the 
system. In other words, there is long-lived entanglement (stationary state entanglement) 
in the presence of the dissipation of cavities, implying that the initial entanglement of two 
atoms could be partially preserved even they are put into dissipative cavities. We find the 
long-lived entanglement depends on the initial state of the two atoms and the ratio of the 
decay rate of the cavities to the atom- field coupling constant. Finally, we discuss the co- 
herence of the two atoms using the linear entropy. Our results show that the linear entropy 
of each atom at any time is equal to the initial linear entropy, that is, the coherence of 
each atom is preserved. We find that the initial coherence (entanglement) of the two atoms 
could be preserved (partially preserved) in the dispersive limit even if they are put into two 
dissipative cavities. 

The present paper is organized as follows. In Section II, we obtain an explicit analytical 
solution of one atom interacting with a dissipative cavity in the dispersive limit. In section 
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III, we consider a quantum system consisting of two atoms within two spatially separated 
cavities. In Section IV, the entanglement dynamics and coherence of the two two-level atoms 
are investigated by employing the concurrence and linear entropy, respectively. Finally, we 
summarize our results in Section V. 



II. ONE ATOM IN A DISSIPATIVE CAVITY 



We first consider a quantum system consisting of a two-level atom interacting with a 



single-mode cavity. Under the electric dipo 
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e and rotating wave approximation, the Hamil- 
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tonian of the present system is {h = 1) 

H = ua'a+'^az + g{a^cr_ + aa-i-), (1) 

where g is the atom-field coupling constant, = |e)(e| — \g){g\ and a± are the atomic spin 
flip operators characterizing the effective two-level atom with frequency ujq. The symbols 
|e) and \g) refer to the excited and ground states for the two-level atom. Here, and a 
are the creation and annihilation operators of the field with frequency u, respectively. The 



dispersive limit is obtained when the condition |A| = Icuq — tul ^ y/n+Tg is satisfied for 
any relevant n. Then, the interaction Hamiltonian g{a'a_ +aa+) can be regarded as a small 
perturbation and the effective Hamiltonian of the model can be recast as 25 1 



He = uja)a + —a, + VL[{a)a + l)|e)(e| - a)a\g){g\], (2) 
with Vt = g'^/A. In the interaction picture, the interaction Hamiltonian is 

V = Q[{a^a + l)\e){e\-a^a\g){g\]. (3) 

We assume the two-level atom interacting with a coherent field in a dissipative environment. 
This interaction causes the losses in the cavity which is presented by the superoperator 
T> = k{2a ■ — a'' a ■ — ■ a'' a), where k is the decay constant. For the sake of simplicity, we 
confine our consideration in the case of zero temperature cavity. Then, the master equation 
that governs the dynamics of the system can be written as follows 

= —i[V, p'] + k{2ap'a) — a^ap' — p'a)a) 

= -i[V,p'] + k{2Mp' -Up' - Cp'), (4) 



where p' is the density matrix of the atom-field system and A4, TZ, and JC are defined by 

A4p' = (a • a^)p' = ap'a\TZp' = {a)a-)p' = a^ap',Cp' = {■a^a)p' = p'a^a. (5) 

Here the superoperators a-, -a, a^-, and -a^ represent the action of creation and annihilation 
operators on an operator 

(a-)o = ao, {■a)o — oa, (a^-)o = a^o, {■a^)o — oa^ . (6) 
We can express the density matrix in the following form 

pit) = P'eeit) ® \e){e\+Pgg{t) ® \g){g\ + p',g{t) ® \e){g\ 

+p',M®\9){el (7) 

where p^^'s are defined as p[j — {i\p'\j), p\j = pj^, i,j — e,g. A straightforward calculation 
shows that 

^^'-^^^ ^ {-in{n-/:) + k{2M-n-jC)}p'jt)^/:Ueit), 



dt 

= mn-£) + k{2M-n-£)}p'Jt)^£,,p'Jt), 

= {-in{n + c + i) + k{2M-n-c)}p',g{t) 

P'geit) = Pl,{t). (8) 

It is easy to check that the superoperators Jv[, TZ and C satisfy the relations 

[n,M\ = [i^,M\ = -M,[n,c] = Q. (9) 

It is worth noting that \TZ-\- C,M\ — —2A4, the superoperators TZ + C and M. form a shift 
operator algebra. Thus we have the expansion of the exponential of a linear combination of 
TZ + CemdM 

^jCeet ^ ^{e'^''*-i)M ^-{m+k)fR^{in-k)t/: 

g£egt _ g-img[e2('"+'=)*-l]A</(jQ+fc)g-(in+fc)t7^g-(in+fe)t£ 

Suppose the atom is initially prepared in the state j 1 g^j^d the field is initially 

prepared in the coherent state jo;) = c"'"'^/^ ^ifl^) with a being a complex number. 
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Here, \n) is the Fock state with a)a\n) = n\n). Therefor the initial state of the atom-cavity 
system is 

p'(0) = ( ^\a){a\ 

= {Ue){e\+^b\9){9\ + ^c\e){g\ 

+CI^)(e|)®|«)(«|. (11) 
Comparing Eq.© with Eq. ffTT]) . one get 

p'ee(O) = ^a\a){a\,pgg{0) = ^b\a){a\, 

Pe,(0) = ^c\a){a\,p'gM=pliO)- (12) 

Combing Eqs.(lO), f|T2l) and the relation p'^jit) = e^'^Vjj(O)) "we find that the matrix elements 
Pijit) at time t is given by 

P'eeit) = ^a\a+it)){a+{t)\, 

P'ggit) = ib\a^{t)){a.{t)l 
P'egit) = Uma^{t)){a.{t)lp'g,{t) = p}g{t) 
|a±(t)) = |ae-(™)*), 

fit) = exp {-iQt + \a\\e^^''' - 1)} 

xexp{J^[l-e-('=+^^)1}. (13) 

The reduced density matrix of the atom is obtained by tracing out the variables of the field 
from the above density matrix 

P'atomit) = '^a|e)(e| 

+ [Um{t)\e){9\ + h.c], 
m = (14) 

where h.c denotes the Hermitian conjugate. 

III. TWO ATOMS WITHIN TWO SPATIALLY SEPARATED CAVITIES 



In this section, we consider a quantum system consisting of two noninteracting atoms 
each locally interacts with its own coherent field of a dissipative cavity. The two atoms are 
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initially prepared in entangled states and the two fields are prepared in coherent states 
and \a2)- For the sake of simplicity, we assume ai = a2, the decay rates of the two cavities 
are equal, and the atom-field coupling constants are the same. The schematic picture of the 
present model is presented in Fig.l. 



A. Reduced density matrix of two atoms 

Using the method introduced by Bellomo, France, and Compagno [isl, we can obtain the 
reduced matrix density of two atoms conveniently. In the basis |1) = |ee), |2) = \eg), |3) = 
\ge), |4) = \gg), and using Eq. IHM . the reduced density matrix p(t) for the two-atom system 
is 

Pii{t) = pii{0), 
Puit) = /(t)A(t)pi2(0), 
PisW = /(t)A(t)pi3(0), 
PmW = [/(t)A(t)]V4(0), 

P23(t) = |/(t)A(t)|V23(0), 
P24(t) = /(t)A(t)p24(0), 

P34(t) = /(t)A(t)p34(0), (15) 

with Pij{t) = p*i(t) and i = 1,2,3,4. We would like to point out that the above procedure 
allows us to obtain the reduced density matrix of the two-atom system for any initial state 
of the two atoms. 



B. Extended Werner-like states 



states 


Q defined by 




P* = 






P* = 






1$) = 


p\eg) + jy\ge), 




|vl>) = 


p|ee) + i^\gg), 





(16) 
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where p is a real number which indicates the purity of initial states, J is a 4 x 4 identity 
matrix, fi and u are complex numbers with + |z/p = 1. The parameter p is 1 for pure 
sates and for completely mixed states. It is worth noting that the EWL states belong to 
the class of the 'X' states. Explicitly, if the density matrix of a quantum state is of the form 
/ pii pi4 ^ 

P22 P23 

Pis P33 , (17) 

Pm P44 

V / 

then it belongs to the class of the X states. 

The EWL states have the following advantages. First, we can easily find that if the initial 
state of two atoms is X structure, then the reduced density matrix at arbitrary time t is still 
X structure under the single atom evolution determined by the master equation (4) (in the 

basis |1) = I + +), |2) = I + -), |3) = | - +), |4) = | )). Second, the EWL states allow 

us to clearly show the influence of the purity and the amount of entanglement of the initial 
states on the entanglement dynamics of two atoms simultaneously. The purity of the EWL 
states are dependent on the parameter p and the amount of the entanglement of the EWL 
states are related to p and u. If p = 1 the EWL states reduce to the Bell-like states |$) and 
In the case oi p = 1, p = z/ = 1/V2 the EWL states become the Bell states while in the 
case of p = they are the maximally mixed states. Third, as we will see, it is easy for us 
to calculate the entanglement dynamics of two atoms which are initially prepared in EWL 
states. 

IV. ENTANGLEMENT AND COHERENCE 

In this section, we will analyze the entanglement dynamics and coherence of two atoms 
by employing the concurrence and linear entropy, respectively. In order to study the entan- 
glement of above system described by density matrix p, we adopt the measure concurrence 



which is deflned by 26 1 



C = max{0, Ai - A2 - A3 - A4}, (18) 

where the A, (i = 1, 2, 3, 4) are the square roots of the eigenvalues in decreasing order of the 
magnitude of the "spin-flipped" density matrix operator R = p{ay ® ay)p*{ay ® ay) and ay 
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is the Pauli Y matrix, i.e., cr„ = ( . Particularly, for X structure states defined by 

V V 

Eq. (17), it is easy to find the concurrence is 

C{t) = 2max{0, jpssWI - Vpn{t)pu{t), \pu{t)\ - Vp^MpM) (19) 
Combing Eqs. (15), (16), and (19), we find the concurrence of the two atoms, which are 
initially prepared in EWL states (p$ or p$), are 

C^{t) = 2max{0,|/(i)A(i)|V23(0)|-A/M0W0)}, 
C^{t) = 2max{0,|/(i)A(i)|>i4(0)|- V^^^MMO)}- 

(20) 

Inserting the initial state of two atoms of Eq.(16) into the above equations leads to 

C{t) = C^{t) = C^{t) 

= max{0,2p|/(i)A(i)|>i.| - (21) 

We find the concurrence C$(t) and Cq,{t) are the same if the parameter /i and u of p$ and 
are equal, respectively. It is clear that the entanglement of the two atoms depends on 
the initial state of the them, i.e., C{t) relies on the parameters p, /i, and u. 

In Fig. 2 and Fig. 3, the concurrence C{t) is plotted as a function of the dimensionless 
scaled time D,t and the parameter fc/Q for different values of p. From these two figures, one 
can easily find that the entanglement of the two atoms suddenly goes to zero and stays zero 
for a finite time, i.e., the ESD phenomenon appears in the present system. However, the 
entanglement cannot revive completely due to the dissipation of the cavity fields. The lower 
panels of Fig. 2 and Fig. 3 are the contour plots of the concurrence, where the severe shading 
areas indicating the two atoms are completely disentangled. Comparing the upper panels 
of Fig.2 and Fig. 3, one can clearly see that the entanglement dynamics relies heavily on the 
purity of the initial state which is represented by the parameter p. For instance, the maximal 
value of concurrence in Fig.2 is about 0.7 while the maximal value of concurrence in Fig. 3 
is only about 0.4. The contour plots of these figures clearly show that the areas showing 
entanglement significantly decrease with the decrease of the purity of the initial state. To 
show this more intuitively, we plot the concurrence as a function of the dimensionless scaled 
time fit for several values of p in Fig.4. Also, the entanglement dynamics depends on the 
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initial state of the field \a) as one can see from the upper panel of Fig. 5. The concurrence 
decreases with the increase of the amplitude of a. For instance, in the case of k/fl — 0.01, 



a = 2, the ESD phenomenon could appear. 

In order to see the time limits (long-term) behavior of the entanglement we plot the 
concurrence as a function of the dimensionless scaled time fit {0 ^ fit ^ 500) for a — 0.5 
(solid line), a — 1 (dashed fine), and a — 2 (dotted line) in the lower panel of Fig. 5. From 
this figure, we see that the entanglement of the two atoms survives for a long time if the 
mean photon number of the fields jap has small values (see the solid and dashed lines of 
Fig. 5). The entanglement of the stationary state p(oo) depends on the initial state of the two 
atoms and the ratio of the decay rate of two cavities to the atom-field coupling constant. If 
the initial state is completely mixed (p = 0), that is, the two atoms are disentangled initially, 
then no entanglement will be generated in the system. 

In Fig. 6, we have presented the time-evolution of entanglement as a function of kt and 
fl/k. One can see clearly that the entanglement decreases gradually as kt increases for a 
given time. In principle, by using the explicit expression of the concurrence, one can find 
exactly the threshold value after which a long-lived entanglement can be obtained. It is 
easy to see that the oscillations of the concurrence gradually disappears as the time kt 
increases. Actually, from Eq.(21) one can see that C{t) — c (c is a constant) leading to 
the following condition: kt > ktc which can be obtained from the following equation(we set 



then, the entanglement has no chance to be oscillating when kt exceeds a threshold value 
ktc (from Fig. 6, ktc ~ 7). Here we should point out that the effective threshold value of the 
decay time corresponding to the transition of the entanglement from oscillating regime to 
fully constant vahic is independent of the interaction time. From the above discussions, one 
can observe that the initial entanglement of the two atoms could be partially preserved in the 
present system which is useful for quantum information processing and quantum memory. 

In order to show the dependence of the long time behavior of the entanglement of the 
atoms on the parameter a and the purity of the initial state p, we plot the concurrence as a 
function of a and p with t — 500 in Fig. 7. From Fig.7, we observe that the stationary state 



jJi = V = 



l/\/2, p = 0.9, and a = 0.5, there is no ESD. However, in the case of a = 1 or 



li^v^ 1/V2) 




(22) 
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entanglement of the atoms decreases with the increase of the parameter a while increases 
with the purity of the initial state p. This observation is consistent with the previous 
discussions. We have also carefully investigated the entanglement dynamics of the system 
if only one cavity is subjected to decay radiation field and find that the entanglement will 
also be partially preserved. 

Now, we investigate the coherence of the atoms by employing the linear entropy defined 



by 



28| 



Sit) = l-Tr[p\t)]. (23) 

The linear entropy 5* is zero for pure states and 1 for completely mixed states. Tracing out 
the degrees of the freedom of the second atom from Eq. (17) and using Eq. (15), we get the 
reduced density matrix of atom 1 

Pi{t) = [pii(0)+P22(0)]|e)(e| 

+ [P33(0)+P44(0)]|^)(^|. (24) 

Thus the linear entropy of atom 1 is 

Si{t) = Si{0) = 1 - Tr[pl{t)] = 2|/iz/p > 0. (25) 
Similarly, for atom 2, we have 

S2{t) = ^2(0) = 1 - Tr[pl{t)] = 2\pu\^ > 0. (26) 

From the above equations, we find that the coherence of the atoms is preserved in the present 
model, that is, the subspace of the atoms are decoherence-free subspaces 



V. CONCLUSIONS 



In the present paper, we have investigated the entanglement dynamics and coherence of 
a quantum system formed by two two-level atoms interacting with two spatially separated 
and dissipative cavities in the dispersive limit. With the help of superoperator method, 
we obtained an explicit expression of the reduced density matrix of the two atoms and 
calculated the entanglement and coherence of them by employing the concurrence and linear 
entropy, respectively. In a short-term regime, the ESD phenomenon could be obtained and 
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the entanglement of the two atoms decreases with the time development. However, in the 
long-term, the entanglement of the two atoms tends to a fixed value. This value depends 
only on the initial state of the two atoms, the atom-field coupling constant, and the decay 
rate of fields. Particularly, we found that there is long-lived entanglement (or stationary 
state entanglement) in the presence of the dissipation of the fields. In other words, in the 
dispersive limit, the initial entanglement of the two atoms within two dissipative cavities 
could be partially preserved for a long time. 

Finally, we calculated the coherence of the two atoms using the linear entropy. Our results 
show that the coherence of each atom can be preserved. Thus the subspace of each atom is a 
decoherence-free subspace. In a word, the initial coherence (entanglement) of the two atoms 
could be preserved (partially preserved) in the dispersive limit even if they are put into two 
dissipative cavities. This feature is useful for quantum information processing and quantum 
memory. It is interesting to extend the present system to a many-qubit system. Note that 
Eq.(15) is not suitable for the many-qubit system. However, the method introduced by 
Bellomo, France, and Compagno 13|] is still applicable. Thus, one can investigate whether 
the multipartite entanglement of the many-qubit system can be preserved. 
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dissipative cavity dissipative cavity 



FIG. 1: The schematic picture of two two-level atom within two spatially separated and dissipative 
cavities. The two atoms are initially prepared in EWL states. Each atom is put into a single-mode 
cavity and interacts with its own cavity field locally. Note that there is no direct interactions 
between two atoms once they have been put into cavities. 




FIG. 2: The concurrence of two atoms is plotted as a function of the dimensionless scaled time fii 
and parameter k/Q. with a = 1, ^ = u = l/\/2) and p = 0.8. 




FIG. 3: The concurrence of two atoms is plotted as a function of the dimensionless scaled time fii 
and parameter k/Q. with a = 1, fi = v = 1/^/2, and p = 0.6. 
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FIG. 4: The concurrence of two atoms is plotted as a function of the dimensionless scaled time 
Clt with k/Q = 0.01, a = 0.5, ji = v = l/-\/2 for p = 0.8(solid line), p = 0.6(dashed line), and 
p = 0.5(dotted line). Upper panel: The short-term behavior of the entanglement between two 
atoms. Lower panel: The long-term behavior of the entanglement between two atoms. 
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at 

FIG. 5: The concurrence of two atoms is plotted as a function of the dimensionless scaled time 
fit with k/Q = 0.01, p = 0.9, jj, = v = l/\/2 for a = 0.5(solid line), a = l(dashed line), and 
a = 2 (dotted line). Upper panel: The short-term behavior of the entanglement between two 
atoms. Lower panel: The long-term behavior of the entanglement between two atoms. Comparing 
the upper panel with the lower panel, one can observe that, in the short-term, there is ESD and the 
maximal values of entanglement decrease with time. However, in the long-term, there is long-lived 
entanglement (see the solid and dashed lines of the lower panel). The entanglement decreases with 
the increase of the mean photon number of the field represented by lap. If the parameter |ap is 
larger enough, there is no long-lived entanglement (see the dotted line of the lower panel). 




FIG. 6: The concurrence of two atoms is plotted as a function of the dimensionless scaled time kt 
and parameter ^/k with a = 1, fj, = v = l/\/2 for p = 0.8. 
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FIG. 7: The concurrence of two atoms is plotted as a function of the parameter a and p with 
n = 1, k = 0.01, fi = u = 1/V2, and t = 500. This fi gure clearly shows the dependence of the long 
time behavior of the entanglement of the atoms on the parameter a and the purity of the initial 
state p. From this figure, we see that the stationary state entanglement of the atoms decreases 
with the parameter a while increases with the purity of the initial state p. 



